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The Stellar-Disk Electric (Short) Circuit: 
Observational Predictions for a YSO Jet Flow 
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Abstract 

We discuss the star-disk electric circuit for a young 
stellar object (YSO) and calculate the expected torques 
on the star and the disk. We obtain the same disk mag- 
netic field and star-disk torques as given by standard 
magnetohydrodynamic (MHD) analysis. We show how 
a short circuit in the star-disk electric circuit may pro- 
duce a magnetically-driven jet flow from the inner edge 
of a disk surrounding a young star. 

An unsteady bipolar jet flow is produced that flows 
perpendicular to the disk plane. Jet speeds of order 
hundreds of kilometres per second are possible, while 
the outflow mass loss rate is proportional to the mass 
accretion rate and is a function of the disk inner radius 
relative to the disk co-rotation radius. 



Keywords Accretion disks 
Outflows 



Stellar magnetospheres 



1 Star-Disk Magnetic Interaction 

Most studies of the electro-magnetic interaction be- 
tween young stars and their nascent discs undertake 
their analysis using the standard magnetohydrody- 
namic (MHD) approximations (Uzdensky 2004). This 
approach has the advantage of describing the Lorentz 
force in terms of magnetic fields and allows an analysis 
that can ignore electric fields and currents. In this pa- 
per, we examine the star-disk electric circuit to see if 
we not only obtain the same answers as standard MHD 
analysis, but also if we can gain new insight into how 
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Fig. 1. — Current flows in the inner section of a 
star /disk circuit, ju - disk current, jm - field aligned 
stellar magnetosphere current, jc - field aligned coro- 
nal currents Rt - inner disk truncation radius, R co - 
co-rotation radius. If the stellar magnetic field pointed 
in the opposite direction, to that shown in the figure, 
then the direction of the current flows would reverse. 

accretion disks and bipolar jets may be related in young 
stellar systems. 

It is assumed that a star is rotating with an angular 
frequency f2*, it has a dipole magnetic field B*, and 
that the magnetic moment /x* is aligned with the ro- 
tational axis (/*„ = f2*), which is perpendicular to the 
plane of the disk (Fig. [T]). The direction of B* is such 
that the z component is negative as it passes through 
the accretion disk. 
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The stellar magnetic field truncates the disk at a 
radial distance R t from the centre of the star, where 
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with /io the permeability of free space, M a the disk mass 
accretion rate, is the magnetic field strength at 

the surface of the star, i?* is the radius of the star, 
M* the stellar mass and G the universal gravitational 
constant. 

The co-rotation distance, i? co , is the radial distance 
from the star where the angular frequency of the stel- 
lar magnetic field (ss f2«) equals the Keplerian angular 
frequency of the disk Qk (r): 
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with P* the rotational period of the star. f2x(f ) is given 
by the equation 
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The relative difference in angular velocity between 
the disk and the co-rotating stellar magnetic field gen- 
erates a current within the disk. In Fig. [TJ we show 
a section of the stellar/disk circuit. Here the current 
density generated within the disk, ju, travels along the 
inner stellar magnetic field lines, jm, and then returns 
to the disk via the outer stellar magnetic field lines in 
the corona above the disk, jc- Only the current flow 
between R t and R co is shown. The full star-disk circuit 
is shown in Bardou & Heyvaerts (1996). 



2 The Star-Disk Electric Circuit 

In Fig. [21 we schematically depict the current flows in 
and around the inner region of the disk. The electric 
field in the disk, Erj(r), is given by (Liffman and Bar- 
dou 1999) 



E n (r)=r(Sl K (r)-n,)B„(r)(' ) 



with \B. t \(r) r* B.(R.) ( — 
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where, from Fig. [IJ B rz (r) — — \B*g(r)\. 

For a disk with finite conductivity, up, the induced 
electric field drives a radial current in the disk with a 
current density of the form 



JdM = -<TD(r)rtt K (r) 
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This radial disk current generates a toroidal magnetic 
field in the disk. The equation for which is (Campbell 
1992): 



B$ = ftQcrrzfo* - Vl K {r))B z {r) . 



(7) 



Campbell used standard MHD to derive Eqn j7]), but 
the same result is also obtained from the current flow 
model of Fig. H (Liffman & Bardou 1999). 
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Fig. 2. — Current flows and magnetic fields near or 
in the disk. The poloidal field, B z , is the section of 
the stellar magnetosphere that interacts with the disk. 
The toroidal disk field, Bp^, and the current flows are 
generated by the relative motion between the disk and 
B z . In this case, jr> is the disk current density, ju 
is the magnetospheric current density that travels be- 
tween the inner edge of the disk and the star, while j z 
is the z component of jc - the current between the star 
and the disk. 

To compute the field aligned current, j z i.e., the z 
component of jc that enters the disk (Fig. \2§ , we apply 
the steady state form of conservation of electric charge 
(V- j = 0), which implies 



Jz = 



-1 
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where (3 is a non-dimensional parameter with the defi- 



nition (Matt and Pudritz 2005) 
/3 _1 = ^a D rhVl K (r) — 2.8 x 

(io- 7 Sin 1 



(9) 
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with h being the scale height of the disk. 

We denote by Im the total current from the top or 
bottom half of the inner disk (i.e., the inner edge lo- 
cated at the truncation radius, R t ) that travels along 
the stellar field lines to the star (the corresponding cur- 
rent density, jm is shown in Fig. [I]), the magnitude of 
Im is given by (using Eqn ©) 
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A representative value for the magnitude of Im is given 
by 
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3 Disk-Star Torque 




Substituting Eqns ((6]) and j9|) into Eqn (TT2"]) gives the 
gradient of the torque exerted by the stellar magnetic 
field onto the disk: 
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Substituting Eqns © and © into Eqn (O gives 



Ar 



Mo 



(14) 



The same equation has been derived via standard MHD 
analysis (Clarke et al. 1995). This suggests that the 
equations for the currents and magnetic fields, as given 
here, have the correct form. 



4 The Short Circuit Model 

We now assume that a portion of the inner field-aligned 
current, jM, (Fig. H]) short circuits and produces a ra- 
dial current. A quantitative discussion of such trans- 
field current flows is given in Chapters 4 and 7 of Brekke 
(1997), where it is shown that transfield currents, such 
as gravitational drift currents, regularly occur in the 
Earth's ionosphere and magnetosphere. 
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Fig. 3. — The interaction of an assumed bipolar stellar 
magnetic field, B and the disk produces current flows, 
J , (denoted by fat arrows) that, in turn, create J x B 
forces which act upon the star and the disk. For the sys- 
tem shown, the general rotation is anticlockwise when 
viewed from above. The Q and (^) symbols represent 
the Lorentz force pointing in the direction towards and 
away from the observer, respectively. 

To determine the torque(s) on the disk, we consider 
an annulus of the disk, which has a radius of r, thickness 
Ar and height 2h. The volume, AT^, of the annulus is 
4iTrhAr and it feels a torque 

Ar = |r x (j D x B* z )\ AV = Ai:r 2 hj D B, z Ar (12) 



Fig. 4. — The z component of the stellar gravita- 
tional field, g z , interacts with the toroidal magnetic 
field above the disk to produce a radial drift current 
with a velocity Vorift ■ 

To illustrate how, for example, a radial, gravitational 
drift current could arise from the field configuration 
shown in Fig. we note, from Brekke (1997), that the 
drift velocity, Vd, of a charge, q, subject to a force, F±, 
perpendicular to a magnetic field, B, is given by 



F±x B 

qB 2 



(15) 



Using Eqn (flBf we schematically show, in Fig. [H 
how the cross product of the "wound up" toroidal field 
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(i.e., Bm4>, produced by the disk/star current flow, 
Im) and the z component of the stellar gravitational 
force can produce a radial, gravitational drift of pos- 
itively charged particles (and hence a current) above 
and below the disk. Other particle drifts are also pos- 
sible, the gravitational drift case is shown for purposes 
of illustration. 

It is presumed that this hypothetical 'short-circuit' 
region has an inner radius of r = R t and an outer ra- 
dius, i? , where R < R co . The radial transficld cur- 
rent, (j r ), can interact with the toroidal field, Bm</>- 
The subsequent j r x Bm4> Lorentz force is in the cor- 
rect direction to power a flow away from the disk. 
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Fig. 5. — Even in a highly-conductive plasma, trans- 
field currents can flow between the stellar magnetic field 
lines above and below the accretion disk between R t 
and R . The interaction between the transfield currents 
and the toroidal fields give rise to the j x B Lorentz 
forces, which drive the outflow 

For the jet flow shown in Fig. 03 the flow can only 
escape the stellar magnetosphere when: 
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where, p is the gas mass density, v the wind speed and 
it is assumed that the main part of the flow occurs at 
the inner edge of the disk (r — Rt). To find the required 
values of p and v, we note that the velocity of the flow 
has to be of order the escape speed: 



r y (r/0.05AU) 

Combining Eqns (fTH)) and (fTT|) gives 

BJr) 2 r 
p ~y^r- = 1-86 x 1(T 13 x 



(g*CR,)/0.01T) 2 (i?,/2R e ) 6 
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The existence of such a 'break-out' energy suggests the 
possibility of a pulsatile jet flow. 

4.1 Jet Exhaust Speed 

Applying Amperes Law to a thin slice (thickness dz) of 
a region above the disk gives 



dB 4 
dz 



(r,z) = -fj, jr(r,z) 



(19) 



As an illustrative example, we will assume a constant 
radial current in z for z € [zo, Zt], where zq and zt are 
the bottom and top, respectively, of the outflow acceler- 
ation region. We denote this constant, z independent, 
radial current density by J r . By also assuming that 
Bc/)(r, zt) = 0, we can solve Eqn (fT9"| to obtain 



B$(r, z) = p (z T - z )J r (r) 1 



Z - Z 



(20) 



zt - zo , 

The total radial current, I r (r), is given by 
I r (r) = 2-,rr(z T - z )J r {r) . (21) 
Combining Eqns l[2U]) and (PT]l gives 

B<j,{r, z) = — 1 , z € [z , z T \ . (22) 
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Liffman & Siora (1997) obtained a Bernoulli equa- 
tion for the flow: 



7 \ p B 2 GM* GM*r 



7 - 1 / p p, p ^/r 2 + z 2 2r 2 



■■£, (23) 



with £ the constant specific energy of the streamline 
flow, 7 - the ratio of specific heats, B is the toroidal 
magnetic field, p - pressure, tq - the initial value of r, 
and v — yfv 2 + v 2 . 

Using Eqns (|2"2"|) and (|23|l one can obtain an expres- 
sion for the exhaust speed of the jet flow, v e : 
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(24) 
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where po is the gas density at the base of the flow. 
The representative values for po, I r and ro, as given in 
Eqn pi]) , are obtained from Eqns {IS]). (fTT|) . and ((TJ, 
respectively. 
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4.2 Radial Size of the Outflow Region 

We denote by R the distance from the star, where 
the integrated current density entering the top half of 
the disk (j z - as depicted in Fig. [5]) is equal to the 
integrated current density returning to the disk via the 
stellar magnetosphere (j'm and j r as in Fig. 

Let I z (r) denote the total vertical current entering 
the top half of the disk from the inner edge of the ac- 
cretion disk, to a distance r from the star: 



h{r) = / j z {r)2-Krdr 
JR t 



(25) 



Substituting Eqns © and {5} into Eqn ([2"5]) implies 
4(r) = - 27riWJ?t) x 



(26) 
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If we consider the case where all of the Im (Eqn (jTOJ) ) 
current short circuits via the radial transfield current j r 
and then back to the disk via the field-aligned current 
j z then, for this case, R is determined by equating Im 
and I z {f), and specifying that r — R . This condition 
implies that 
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We can numerically solve for R a in the above equa- 
tion and obtain values for (R a — R t )/R co , which is the 
normalized length of the outflow region in the inner ac- 
cretion disk. These results are shown in Fig. where 
it can be seen that as the inner disk approaches the 
star (R t — > i?*) the width of the outflow acceleration 
region decreases (R — ► Rt). Similarly, as the inner 
disk approaches the co-rotation radius the width of the 
outflow acceleration region decreases to zero. The max- 
imum width of the acceleration region occurs when the 
inner disk radius is approximately half that of the co- 
rotation radius. This behaviour is represented schemat- 
ically in Fig. [HI , where we note that the outer radius of 
the outflow acceleration region is always less than the 
co- rotation radius (R a < R co )- 

4.3 Mass Ejection Rate 

From the conservation of mass, the mass ejection rate 
of an outflow, M w , is 
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Fig. 6. — The length of the outflow active region of the 
inner disk (R —R t ) as a function of the inner truncation 
radius, Rt, of the disk, where both quantities are nor- 
malized to the co-rotation radius, R co . The widths of 
the outflow at points A, B and C are depicted schemat- 
ically in Fig. [H 
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Fig. 7. — The ratio of outflow mass rate, M w , to the 
mass accretion rate onto a star, M a , versus the ratio 
of the inner disk truncation radius to the co-rotation 
radius (Rt/R co ). The mass flow rates of the outflow at 
points A, B and C are depicted schematically by the 
length of the arrows in Fig. [8] 

where p, v and A are, respectively, the density, speed 
and cross-sectional area of the outflow. Noting that 



A = n(R 2 - R 2 t ) 

and using Eqns {TTJ) , (HH and HJ, Eqn 
form 
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(29) 
has the 
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Using Eqns (|3"D|) and (f2"7) . we can compute the ratio 
M w /M a as a function of R t /R co . These results are 
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Fig. 8. — A schematic depiction of the mass outflow 
rate and the radial size of the outflow acceleration re- 
gion as a function of the inner disk truncation radius, 
Rt, and the co-rotation radius, R co - indicated by the 
line. The length of the arrow represents the outflow 
mass rate, M w , while the width of the arrow indicates 
the actual, relative size of the outflow acceleration re- 
gion. Case A: R t is small relative to R co and M w ~ M a . 
B: R t ~ R co /2, the outflow acceleration region is at its 
broadest and M w ~ M a /2. C: R t ~ R co and M w — ► 0. 

shown in Fig. [7J where it can be seen that when R t — 
R co the high-speed outflow shuts down. On the other 
hand, for a fixed co-rotation radius, the mass outflow 
rate increases and approaches the mass accretion rate 
as the inner disk radius approaches the surface of the 
star. This behaviour is represented schematically in 
Fig.H 

The observed values for the mass outflow and accre- 
tion rates are only known to order-of-magnitude values: 
M w /M a ~ 0.1 (Calvet 1997), while the modeling of T 
Tauri observational data gives: Rt/R co ~ 0.6 to 0.8 
(Kenyon et al. 1996). From Fig. [3 the corresponding 
range of values for M w /M a is M w /M a m 0.33 to 0.13 . 
These values are consistent with the observational val- 



short circuit may be due to a gravitational drift cur- 
rent, which might be of sufficient magnitude to gener- 
ate a non-constant bipolar outflow at the inner edge of 
the disk that flows in a direction roughly perpendicular 
to the disk. 

In this model, the irregularity in the outflow arises 
because the outflow has to disrupt the stellar magnetic 
field for it to escape the stellar-disk system. The model 
predicts that the mass outflow rate is proportional to 
the total mass accretion rate in the disk. The mass out- 
flow rate is also dependent on the position of the inner 
edge of the disk relative to the disk co-rotation radius. 
If the position of the inner edge of the disk is equal to 
the disk co-rotation radius, then there is little or no 
outflow. As the inner edge of the disk approaches the 
star (assuming a constant co-rotation radius) then the 
proportion of material going into the outflow increases, 
while the proportion of material accreting onto the star 
decreases. 



References 

Bardou A., Heyvaerts J.: A&A 307, 1009 (1996) 
Brekke A. : Physics of the Upper Polar Atmosphere, Wiley 
(1997) 

Calvet N.: In Reipurth B., Bertout C. (eds) IAU Symp. 182, 
Herbig-Haro Flows and the Birth of Low Mass Stars, 
p. 417. Kluwer, Dordrecht (1997) 

Campbell C.G.: Geophys. Astrophys. Fluid Dynamics 63, 
179 (1992) 

Clarke C.J., Armitage P.J., Smith K.W., Pringle J.E.: MN- 

RAS, 273, 639 (1995) 
Kenyon S. J., Yi I., Hartmann L.: ApJ 462, 439 (1996) 
Liffman K., Siora A.: MNRAS 290, 629 (1997) 
Liffman K., Bardou A.: MNRAS 309, 443 (1999) 
Matt S., Pudritz R. E.: MNRAS 356, 167 (2005) 
Uzdensky D. A.: APSS 292, 573 (2004) 



ues. 



5 Conclusions 

The star-disk electric circuit arises due to the inter- 
action of the stellar magnetosphere and the accretion 
disk. We have shown that the electric circuit model 
and standard MHD analysis give the same expressions 
for the disk magnetic field and the torque between the 

disk and the star. 

We examined the hypothetical case where there is This 2 - C olumn preprint was prepared with the AAS IAIftX macros 
a short circuit in the star-disk circuit. This radial v5.2. 



